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Abstract 

We construct a tractable model to describe the rate at which a knotted polymer is ejected from a spherical capsid via a 
small pore. Knots are too large to fit through the pore and must reptate to the end of the polymer for ejection to occur. 
The reptation of knots is described by symmetric exclusion on the line, with the internal capsid pressure represented 
by an additional biased particle that drives knots to the end of the chain. We compute the exact ejection speed for a 
finite number of knots L and find that it scales as l/L. We establish a mapping to the solvable zero-range process. We 
also construct a continuum theory for many knots that matches the exact discrete theory for large L. 
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1. Introduction 



The two basic steps in bacterial infection are the ini t ial inj ection of the viral DNA into a bacteriophage cap- 
sid dRiemer et all 1 19781: iKindt et all |200U IPurohit et all 120031) and then the ejection of the DNA into a host cell 



dlnamdar et all l2006t ICastelnovo and Evilevitchi I2007IK The processes underlying these two steps are both phe- 
nomenologically rich and incompletely understood. In packaging DNA into a capsid, both the strong repulsion 
from the bending of the DNA and the stress caused by the capsid being smaller than the persistence length have 
to be overcome. The pressure associated with this capsid packaging can be in the range of tens of atmospheres 
( Smith et al. , 2001 ; Gelbart and Knobler , | 2009h. These forces that repel the polymer during packaging are overcome 



by a specifically-designed motor protein dlnamdar et al.L 12006) 



The complementary ejection of the DNA into the host cell through small pores in the capsid is driven by the 
osmotic pressu re exerted by the cap s id and the internal stresses that have been built up in the highly-confined DNA 
polymer chain dlnamdar et al. . 2006; G elbart and Knoblen 120091) . Bec ause of its cen t ral im portance in the life cycle 



of viruses, DNA ejection from a capsid has been intensively studied dlnamdar et all |2006). Most theoretical treat- 



ments have investigated the role of microscopic mechanisms that arise naturally from the combined effects of osmotic 
pressure and bending energy. 

Recent experiments and simulations suggest that viral DNA may become knotted (i. e., a closed chain with a kno t 
cannot be smoothly deformed into a circle) as it is being tightly packed in the capsid dArsuaga et al. . I2002L l2005h . 
These knots are appare ntly large enough to prevent the DNA from being completely ejected being from the capsid 
(Mangeno t et all [20051) . If this steric hindrance is fully operative, then knots would have to unravel for complete 
ejection. Because the driving forces are pushing the chain out of the capsid, the only way for complete ejection to 
occur is for the knot to reptate along the chain and ultimately unravel when the end of the chain is reached. This 
mechanism and its role on polymer ejection has recently been explored by num erical simulation of a bead-spring 
model of a flexible polymer together with a coupling to a background solvent (Matt hews et al. . l2009h . This study 
provides detailed, but nevertheless still qualitative results, for the time dependence of the fraction of the chain that 
remains within the capsid, as well as the dependence of the ejection time on the number and type of knots on the 
polymer. 
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In this work, we propose a coarse-grained model to capture the role of knots o n the ejection process (Fig, [p. In 

Urn 



this model, th e reptation of knots is described by the symmetric exclusion process (Ha rris et al.U l965: Liggett, 1985 



Schu5 l2000l) . while the point at the interface between in the interior and exterior of the capsid undergoes a biased 
motion and also satisfies the exclusion constraint. We may view this biased particle as a "shepherd" that pushes the 
knots ("sheep") to the end of the chain where they can unravel. Once all the knots have unraveled, the polymer can be 
completely ejected. 

In the next section, we present the details of this shepherding model. In Sec. [3] we treat the ejection of a polymer 
chain with a small number of knots. For this discrete system, we write and solve the master equations that describe the 
position of the knots along the chain. This solution gives the exact result that the velocity of a strongly-biased shepherd 
equals 1/L, where L is the number of knots. An unexpected and simplifying feature of the strongly-biased limit is 
that the velocity and diffusion coefficient of the shepherd do no t depen d on microscopic rates. We also establish a 
mapping to the solvable zero range process dEvans and Hannev , 120051) . and we use this mapping to obtain further 
exact results for the diffusion of the shepherd. Finally, we investigate a semi-infinite system with a finite density of 
knots (Sec.|4]i. Here we apply a continuum approach to solve for the knot density profile, from which the displacement 
of the shepherd grows as VAf, with a calculable amplitude A. 



2. Model 

To present our model, it is convenient to employ a reference frame that is fixed along the length of the polymer. 
The chain can be divided into a portion that has been ejected from the capsid and the portion that remains in the 
interior. The basic feature of our modeling is that ejection may be hindered by the existence of knots along the portion 
of the polymer that is still insi de the capsid. We assu me that these knots are sufficiently large that they cannot pass 



through the pore in the capsid (Matthe ws et all 12009). Because of this restriction, some other relaxation mechanism 



is needed to allow for complete ejection of a polymer. 
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Figure 1: Schematic illustration of the shepherding process. The motion of the shepherd (T) is biased to the right, while the knots (•) are represented 
as localized particles that hop symmetrically on the integer interval [0, N]. Upon reaching the end of the chain, the knot has unraveled; this event 
is represented as an absorbing boundary at N. 

One such mechanism is provide d by the possibility that knots are not static but move along the polymer chain 



Une sucn mecnanism is provided by tne possibility mat Knots are not static but move along tne polymer cnain 
by reptation (Matt hews et all 20091) . When a knot reptates to the free end of the polymer, it disappears because the 



entanglement of the knot is released. If the chain contains multiple knots, then we assume that they cannot pass 
through each other; instead they must individually reptate to the end of the chain and unravel in sequence. 

Our main focus is the position of the point that separates the ejected portion of the polymer from the portion inside. 
We model the position of this dividing point by an effective particle (a "shepherd") that undergoes biased diffusion 
along the chain, with rates a and y > a to the left and right, respectively. The bias of the shepherd corresponds to 
the polymer being pushed out of the capsid by osmotic pressure and the release of bending energy. The shepherd 
therefore rectifies the diffusive motion of the knots so that the polymer will be ejected. When the shepherd reaches 
the right end of the polymer, the ejection is complete. In the spirit of minimalism, we assume that this bias is constant 
throughout the ejection process; in reality the bias will decrease with time as the ejection of the polymer relaxes the 
stresses that drives the ejection. 

We assume that the chain initially contains a specified number of knots that are randomly interspersed along the 
chain. The ejected portion is necessarily free of knots, while the complementary portion of the chain remains knotted 
(except possibly near the very end of the ejection process). In the spirit of a coarse-grained model, we ignore the 
topology of the knots and replace them by localized defects. Each defect can move equally likely in both direction 
along the chain due to thermal noise, and hence the reptation of the knot is modeled as a symmetric random walk. 
Different type of knots could be described by different diffusion coefficients, however, we restrict our discussion to a 
single type of knot, where all knots have the same hopping rates. In addition to the random motion of the knots, there 
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is also an exclusion constraint because knots cannot pass through each other or pass by the shepherd (corresponding 
to the disallowed event of a knot passing through the capsid pore). The knots can be viewed as a flock of diffusing and 
mutually repelling excitations that are pushed to the end of the interval by an advancing shepherd. When each knot 
reaches the right end of the chain it simply disappears, corresponding to an absorbing boundary condition. When the 
shepherd reaches the end of the chain (Fig. [TJ, the polymer has been completely ejected. 

Because of its bias, the shepherd rectifies the diffusive motion of the knots. This recti fication of the rmal noise 
underlies many microscopic biological processes, such as the motion of motor protein s ([Howard! 1200 11) . and the 
chaperon assisted translocation of polymers across a membrane (iD'Orsogna et all 120071) . In our model of polymer 
ejection, the rectification is mediated by the exclusion constraint of the knots, leading to an exclusion process in 
which a single particle is subject to a bias that acts indirectly on all other particles. A t wo-sided version of this model 
with a biased particl e caged in between unbiased particles has been studied previously ( Burlatskv et all 1992al 1996; 
Landimetal. , 1998); this model mimics, e.g., the field-driven motion of a charged particle that is immersed in a lattice 



gas of neutral particles. It is worth mentioning that various embellishments of idealized exclusion processes have 
helped to und erstand a variety of biological processes and have raised new theoretical questions about the exclusion 



process itself dReichenbach et all 12000: Istukalin et all 120061 bong et all l2007t lAntal etalll2007h . 



3. Discrete Formulation 

We first study the case of a small number of knots by writing the discrete master equations that describe the 
probability distribution for their positions. These master equations turn out to be soluble, from which we can extract 
the speed of the shepherd in the steady state. 

3.1. Single Knot 

To determine the ejection speed of a polymer that contains a single knot, note that before the shepherd first reaches 
the knot, its speed is simply vo = y - a (Fig. [2] left). Subsequently, the knot and the shepherd stay close to each other 
because of the bias of the shepherd toward the knot. As a result, the speed of their mutual advance is less than vo. To 
obtain this speed, we define P(n) as the probability that the number of empty sites between the shepherd and the first 
knot is n. We call this number of vacancies between the two particles as the "gap". For gap size n>0,n can increase 
to n + 1 with rate 1 + a, either by the knot hopping one step to the right or the shepherd hopping one step to the left 
with rate a, respectively. Similarly, for n > 1, the gap size n decreases to n - 1 with rate 1 + y. When n — 0, the gap 
size can only increase with rate 1 + a. 




time time 

Figure 2: (Left) Time dependence of the positions of the shepherd (lower) and a single knot (upper), when the shepherd hopping rates to the right 
and left are y = 1.1 and a = 1, respectively. (Right) Time dependence of the positions of the shepherd (lower) and 10 knots (upper) on a chain of 
length 1000, when the shepherd hopping rates to the right and left are y = 2 and a = 1, respectively. 

Because the shepherd is driven toward the knot for y > a, the distance distribution between these two entities 
reaches a steady state. In this state, the two transition rates (1 + a)P(ri) and (1 + y)P(n + 1) are equal and give the 
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recursion, for n > 0, (1 + a)P(n) — (1 + y)P(n + 1) for n > 0, with solution 



7+1 \ 1 +7 

Here the overall amplitude is determined by the normalization condition 2„>o P( n ) - !■ 

Since the shepherd can always hop to the left (with rate a), but it can hop to the right only if its right neighbor is 
empty (probability 1 - Pq) , the speed of the shepherd for the case of one knot is 

Vi = y(l-P )-a = y\^-a, (2) 

In the limit of an infinitely biased shepherd (y — > oo), Eq. (0 reduces to V\ — 1 . 

3.2. Several Knots 

For L knots their configuration may be specified by the set of separations n = («i,«2, • • • > «lX where n, is the 
number of vacancies to the left of the z* knot. The stationary probability P(n) for a configuration specified by n 
satisfies the equation 

dP(n) 

— '—L =-( a + y + 2L)P(n) + yP(n x + 1) + aP(n\ - 1) 
dt 

l-i (3) 
+ J] [P(n» - l.njh.1 + 1) + P(n k + l,n i+I - 1)] + P(n L - 1) + P(n L + 1) = 0. 

A-=l 

The scalar arguments of P on the right-hand side indicate that only the components of the separation vector that have 
been changed from those in n. Equation (O is valid for configurations in which the separations between consecutive 
knots are all at least one (n^ > 1 for every k). For configurations where some n# = (i.e., zero separation between 
adjacent knots), the relevant separations cannot decrease. In such cases and independently for each k, the following 
terms are missing (i.e., they should be subtracted) from the right-hand side of © 

-(1 + 7)P(n) + aP(n x - 1) + P(n { - 1,« 2 + 1) for k = 1 

-2P(n) + P(n k -i + \,n k -\) + P(n k -\,n k+l + \) for 2 < k < L - 1 (4) 
-2P(n) + P(n L -i + l,n L - 1) + P(n L - 1) for k = L. 

(Here we formally allow n k to take the value -1 for convenience.) By solving the system of equations (01 and for 
a small number of knots, a simple pattern emerges that suggests that the exact solution factorizes and individual gaps 
are distributed according to exponential distributions 

L 

P(n) = Y\a-z k )4 l . (5) 
k=\ 

(The prefactor in (0 assures that normalization is obeyed: 2 P(n) = 1.) Substituting ansatz (0 into the bulk and 
boundary equations, Eqs. (0 and (0, we obtain 

L-l 



a + 7 + 4 = 7z 1 + + +z L +) , (6) 

Zi zl f-£ z k 

as well as 



Zi(l +y) - a + zi , for k = 1 

2Zi =Z*-i +z i+ i, for 2 < k < L-l (7) 
2zl = Zl-i + 1 . for k - L 

When the conditions (0 and (0 are both satisfied the solution is stationary. Equation (0 represents a linear recursion 
with the unique solution 

aL+ 1 +(y-a)(k- 1) 

Zi = r— : ■ (») 

7L + 1 
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This form for Zk also satisfies the bulk equation ||6}, which can be verified by substitution. 

The speed of the shepherd can be calculated in the same way as in the case of a single knot. The shepherd can 
always jump to the left, but it can jump to right only if its right neighbor is empty, that is Vl = yProb(«i > 1) — a. 
The probability that the first gap has size of at least one is 

oo 

Prob(n 1 >l) = (l- Zl )2^ 1 = Zl . (9) 

n, = l 

Hence the speed of the shepherd is 

y — a 

V L = -j— - . (10) 

JL + 1 



is (1 - Zk) z k k 1 and therefore the average distance is 



Equation (0 shows that the probability distribution for the gap size 1% = 1 +raj between the (k - l) st and k ,h knots 

iverage distance is 



/>,>=/(* l 1 - ,-,;).- " ^ jL+l 1 

f-i k l-Zk y-aL-k+1 

Since gaps are independent, the average total size (i.e., the average distance between the shepherd and the L th knot) is 

L 



Zyu + l 
(h)= L H L (12) 



where Hi = Tii<j<Lj 1 are harmonic numbers. Using the asymptotic Hi InL when L » 1 we see that the average 
gap size scales according to 

<£>^— ^LlnL (13) 
y-a 

when the number of knots is large. Similarly, we use the independence of gaps to compute the variance of the total 
size 

L L 

(?) - (If = J] [(4) - <4 > 2 ] = 2 —^— 2 (14) 



which simplifies [upon using ©J to 



P } -(tfJyJ^l) 2 H v-y±±± HL (15) 

\y — a ) y-a 



where = 2i<;<l j 2 - F° r large number of knots — > 7T 2 /6 and therefore asymptotically 

(f)-{() 2 = {—?— \L 2 +0(L\nL) (16) 



6 \y - a 

Hence the relative fluctuations of the total size diminish, albeit slowly as the inverse logarithm of the total number of 
knots: 

~ — (17) 
({) InL 

In the special case of an infinitely biased shepherd (y = oo), the shepherd is always adjacent to the first knot, so 
that n\ — and the stationary probability (0, ([8]l becomes 

L (k-l ' 



k=2 



P(n) = Y] V " (18) 



In this limit, the speed of the shepherd reduces to Vl — l/L and other characteristics of the system similarly simplify. 
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A useful feature of our model is that it can be mapped onto the zero range process (ZRP) dEvans and Hannev 



2005) on a chain of L sites with open boundaries. In this mapping, each knot corresponds to a site in the ZRP, and the 
number of empty sites to the left of each knot corresponds to the number of particles that occupy site k in the ZRP. 
The particles in the ZRP multiply occupy each site, and only the top particle can hop at a rate that is independent of 
the occupancy of the target site. For example, knot k hopping to the right corresponds to a particle in the ZRP hopping 
from site k + 1 to k. The top particle in this corresponding ZRP can hop: (i) both to the left and to the right at rate 
1 from each bulk site; (ii) out of site k = 1 at rate y and into k = 1 at rate a; (iii) either in or out of site k — L at 
rate 1 . The relevant feature of this mapping is that the negative of the current entering the ZRP from the left exactly 
corresponds to the speed of the shepherd. 



Our shepherd model is a special case of the general ZRP that was solved by lLevine et al.l (12005b and Eqs. (Q, {§), 
and ( fTOb are special cases of their solution. In addition to its average speed, we are also interested in the diffusion 
coefficient of the shepherd. Due to the equivalence of the models, the negative time integrated current — J it) of the 
ZRP corresponds to the position x{t) - -J(t) of the shepherd in our model. We can use results of lHarris et alJ |2005) 
for the fluctuations of the time-integrated current. From this work, the large time asymptotic of the Laplace transform 
o f the position of the shepherd is given explicitly by (the la rge-L limit of the expression below was obtained previously 
in lBodineau and Demdald2004 : lwiiland and rIcz! ()2005l) ): 

m^ l -lo,(e-^)= { l^^. (19) 
f->°° t 1 + yL 

From this equation, the speed of the shepherd ( fTOb can be obtained as Vl = c^xl^o, while the diffusion coefficient is 



_ 1 d 2 x 
Dl ~2M 



a + y 

(20) 



, =0 2(1 + yL) 

Again, the result is particularly neat, Dl = (2L)~', in the limit of an infinitely biased shepherd. 
3.3. Ejection Time 

What is the average ejection time of a polymer with a finite density of knots? For concreteness, consider a polymer 
of length N, with a shepherd at site 0, and L equidistant knots initially at x,- = iN/(L + 1), i = l,...,L, If the knot 
density - L/N is sufficiently low (roughly if pco <K y - a), the biased shepherd reaches each knot close to the 
knot's initial position. The shepherd advances with the localized flock of k already-collected knots with asymptotic 
speed Vk = (y — a)/{yk +1), see Eq. ilOi , until the flock reaches the (k + l) st knot. Consequently, the time for the 
shepherd to reach the end of the polymer is 



N yi 1 _ N yL + 2 



L + l^V k 2 y 



(21) 



Simulation results for this time are plotted in Fig.[3]for a polymer of length N = 100 as a function of the number of 
knots when the bias of the shepherd is weak (a — l,y = 2). Our result (|2TT > gives an excellent fit to the data even for 
this short polymer. 

For a finite density of knots = L/N, we also obtain the asymptotic position of the shepherd from OTT) as 



x— — —— t — > ■\[2t~[p<x > a^0. (22) 

\ yp°° 

Notice the lack of dependence on the hopping rate y of the shepherd in the asymmetric limit of a — > 0. 



When a real virus is ejected from the capsid, the pressure should decrease during this process (iGelbart and Knoblei 



2009 ). This effect can be mo deled as a decreasing bias y for the shepherd. Since the bias is so enormous in the capsid 
(Gel bart and Knoblei , 20091) . even the decreased bias can be consider ed as large bias in our m odel. This argument is 



supported by the fairly uniform rate of virus ejection in simulations of Matthews et al. ( 20091) . The ejection time (fJTJ 



for large bias becomes t = NL/2, that is independent of the rates a and y. Knowing this ejection time experimentally 
would give us an estimate for the hopping rates of the knots, (which we used to rescale time with). Unfortunately, we 
do not know of such experiments. 
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Figure 3: Average ejection time for a polymer of length 100 as a function of the number of equidistant knots L. The hopping rates of the shepherd 
are or = 1, y = 2. The line is the function 100 * (1 + L) as given by Eq. C3}. 



4. Finite Density of Knots 

When the number of knots is large, it is more convenient to study the knot density profile by a continuum theory. 
While a finite knot density is perhaps not of direct relevance to the polymer ejection problem, this limiting situation 
leads to an appealing exclusion process that can be solved in a simple way by appl ying a scaling approach to account 
for the effective moving boundary condition caused by the motion of the shepherd (ICranklll987b . 

Let p(x, t) denote the continuum knot density at position x and time f, and let x*(t) denote the position of the 
shepherd (Fig. Initially the shepherd is at x, — and it hops to the right with rate 1 whenever its right neighbor is 
vacant. In the continuum limit, the knot density p(x, t) satisfies the diffusion equation 



dp 

dt 



c?p 
dx 2 



(23) 



The diffusion coefficient D = 1 since the hoppi ng rate of eac h knot equals 1 and the exclusion plays no role in the 
overall density when the hopping is symmetric (ISchiitz , 2000). 



Pc 



P 



X 



Figure 4: Sketch of the knot density as a function of position along the chain in the continuum limit. 



As the shepherd advances, knots accumulate in front of it and the knot density increases monotonically from a 
value as x — > oo to the value 1 at x*. The shepherd advances only when the adjacent site is empty. In the continuum 
limit, the analog of this advancement rule is 

dx t dp | 
dt dx\x=x, 

That is, the probability that there is a gap in front of the shepherd equals the local knot density gradient. Thus we must 
solve the diffusion equation d23l in the region x*(?) < x < oo subject to the initial condition p(x > 0, 0) = and the 
boundary condition p(x* , t > 0) = 1 . 



(24) 
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This type of moving boundary-value, or Stefan, problem (ICranki Il987t iLangen, 1 1987b can be solved by a scaling 
approach. Let us assume that the density profile approaches the scaling form 



p(x,t)=f(a where t=J~ l - 
Using this scaling form, the derivatives of the density are: 

at x„ ox x* ox A x% 

Using these, the diffusion equation d23T l can be written in the separated form 

• 

X * X *~ 

where A is a separation constant. Equation (l24l i for the motion of the boundary gives 

x m x* = -/'(0) , 

which both fixes the separation constant A = — /'(0) and implies that the boundary point advances in time as 

x, = V2Af. 



(25) 



(26) 



(27) 



In scaled coordinates the differential equation for the density profile now becomes f"(£) = — A(l + £)/'(£). 
Integrating this equation, subject to the boundary condition p(x t , t > 0) = 1 (or /(0) =1), yields 



fig) = 1 -A ^drje- A{ '^ ll2) . 
To eliminate the unknown constant A = -/'(0), we use the fact that /(°°) = poo to obtain the relation 

y-»00 

p ra = 1-A djj e- A( * + i 2/2) 
Jo 

= 1 - y/nA/2 e A/2 erfc( JZ/2). 

that (i mplicitly) determines A as a function of poo ; here erfc is the complementary error function ( Abramowitz and Stegunl 
19721) . Using the asymptotic forms of the error function, we can extract the limiting behaviors of the separation con- 
stant A, from which the position of shepherd is given by 



(28) 



I V 2 /Poo Poo 4 

^-^(l-poo) PooTl. 



(29) 



Notice that the limiting expression for poo — > reproduces Eq. d22l that we obtained from the exact discrete solution. 



5. Summary 

We modeled the ejection of a knotted flexible polymer from a capsid in terms of a hybrid dynamical model 
that consists of a single asymmetric exclusion process (representing the boundary point on the polymer between the 
interior and exterior of the capsid) that interacts with a gas of symmetric exclusion processes (the knots). In the 
absence of any external forces, the knots reptate symmetrically along the chain. The effect of the osmotic pressure is 
modeled as the single asymmetric exclusion process that shepherds the knots to the opposite end of the chain where 
their entanglement is released. Once all the knots have unraveled, the polymer can be ejected. 

The ejection speed Vl can be solved exactly for the L + 1-particle system that consists of the shepherd and L knots 
and the result is remarkably simple for a strongly-biased shepherd: Vl = l/L. The underlying exclusion process 
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can also be solved for a finite density of knots by a continuum description. Here the amount of the chain that is 
ejected at time t is not proportional to t, but rather to yfi. This behavior matches that of the discrete solution in the 
limit of a strongly-biased shepherd. In all of our theoretical modeling, we use constant bias throughout the entire 
ejection process. We argue that the motion remains in the strong-bias regime throughout the ejection process due to 
the extremely large initial pressure in the capsid. As long as the system remains in this strong-bias limit, the motion 
of the shepherd and the knots are not influenced by the decrease in the ejection force and it is appropriate to treat the 
bias as constant. 

One final point is that we have assumed that all knots are identical. It may be worthwhile to extend our model to 
the physical more realistic case where knots are non-identical, so that their diffusivities are random variables. 
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